Conformational Free Energy Landscape Research Articles

Conformational free-energy landscapes of a Na+/Ca2+ exchanger explain its alternating-access mechanism and functional specificity

Secondary-active transporters catalyze the movement of myriad substances across all cellular membranes, typically against opposing concentration gradients, and without consuming any ATP. To do so, these proteins employ an intriguing structural mechanism evolved to be activated only upon recognition or release of the transported species. We examine this self-regulated mechanism using a homolog of the cardiac Na+/Ca2+ exchanger as a model system. Using advanced computer simulations, we map out the complete functional cycle of this transporter, including unknown conformations that we validate against existing experimental data. Calculated free-energy landscapes reveal why this transporter functions as an antiporter rather than a symporter, why it specifically exchanges Na+ and Ca2+, and why the stoichiometry of this exchange is exactly 3:1. We also rationalize why the protein does not exchange H+ for either Ca2+ or Na+, despite being able to bind H+ and its high similarity with H+/Ca2+ exchangers. Interestingly, the nature of this transporter is not explained by its primary structural states, known as inward- and outward-open conformations; instead, the defining factor is the feasibility of conformational intermediates between those states, wherein access pathways leading to the substrate binding sites become simultaneously occluded from both sides of the membrane. This analysis offers a physically coherent, broadly transferable route to understand the emergence of function from structure among secondary-active membrane transporters.

In Silico Study of the Early Stages of Aggregation of β-Sheet Forming Antimicrobial Peptide GL13K.

Antimicrobial peptides (AMPs) are of growing interest as potential candidates that may offer more resilience against antimicrobial resistance than traditional antibiotic agents. In this article, we perform the first in silico study of the synthetic ß sheet-forming AMP GL13K. Through atomistic simulations of single and multi-peptide systems under different conditions, we are able to shine a light on the short timescales of early aggregation. We find that isolated peptide conformations are primarily dictated by sequence rather than charge, whereas changing charge has a significant impact on the conformational free energy landscape of multi-peptide systems. We demonstrate that the loss of charge-charge repulsion is a sufficient minimal model for experimentally observed aggregation. Overall, our work explores the molecular biophysical underpinnings of the first stages of aggregation of a unique AMP, laying necessary groundwork for its further development as an antibiotic candidate.

Energy landscapes of homopolymeric RNAs revealed by deep unsupervised learning

Conformational dynamics of RNA plays important roles in a variety of cellular functions such as transcriptional regulation, catalysis, scaffolding, and sensing. Recently, RNAs with low-complexity sequences have been shown to phase separate and form condensate phases similar to lowcomplexity protein domains. The affinity for phase separation and the material characteristics of RNA condensates are strongly dependent on sequence composition and patterning. We hypothesize that differences in the affinities for RNA phase separation can be uncovered by studying sequence-dependent conformational dynamics of single RNA chains. To this end, we have employed atomistic simulations and deep dimensionality reduction techniques to map temperature-dependent conformational free energy landscapes for 20 base-long homopolymeric RNA sequences: poly(U), poly(G), poly(C), and poly(A). The energy landscapes of homopolymeric RNAs reveal a plethora of metastable states with qualitatively different populations stemming from differences in base chemistry. Through detailed analysis of base, phosphate, and sugar interactions, we show that experimentally observed temperature-driven shifts in metastable state populations align with experiments on RNA phase transitions. Specifically, we find that the thermodynamics of unfolding of homopolymeric RNA follows the poly(G) > poly(A) > poly(C) > poly(U) order of stability, mirroring the propensity of RNA to form condensates. To conclude, this work shows that at least for homopolymeric RNA sequences the single-chain conformational dynamics contains sufficient information for predicting and quantifying condensate forming affinities of RNAs. Thus, we anticipate that atomically detailed studies of temeprature -dependent energy landscapes of RNAs will be a useful guide for understanding the propensity of various RNA molecules to form condensates.

In Vitro and In Silico Explorations of the Protein Conformational Changes of Corynebacterium glutamicum MshA, a Model Retaining GT-B Glycosyltransferase.

MshA is a GT-B glycosyltransferase catalyzing the first step in the biosynthesis of mycothiol. While many GT-B enzymes undergo an open-to-closed transition, MshA is unique because its 97° rotation is beyond the usual range of 10-25°. Molecular dynamics (MD) simulations were carried out for MshA in both ligand bound and unbound states to investigate the effect of ligand binding on localized protein dynamics and its conformational free energy landscape. Simulations showed that both the unliganded "opened" and liganded "closed" forms of the enzyme sample a wide degree of dihedral angles and interdomain distances with relatively low overlapping populations. Calculation of the free energy surface using replica exchange MD for the apo "opened" and an artificial generated apo "closed" structure revealed overlaps in the geometries sampled, allowing calculation of a barrier of 2 kcal/mol for the open-to-closed transition in the absence of ligands. MD simulations of fully liganded MshA revealed a smaller sampling of the dihedral angles. The localized protein fluctuation changes suggest that UDP-GlcNAc binding activates the motions of loops in the 1-l-myo-inositol-1-phosphate (I1P)-binding site despite little change in the interactions with UDP-GlcNAc. Circular dichroism, intrinsic fluorescence spectroscopy, and mutagenesis studies were used to confirm the ligand-induced structural changes in MshA. The results support a proposed mechanism where UDP-GlcNAc binds with rigid interactions to the C-terminal domain of MshA and activates flexible loops in the N-terminal domain for binding and positioning of I1P. This model can be used for future structure-based drug development of inhibitors of the mycothiol biosynthetic pathway.

Conformational study into N-alkyl-N′-aryl ureas to inform drug discovery

Ureas are an important functional group in small molecule drugs as well as having wider applications in organic chemistry. Understanding of their conformation is of critical importance for rational design of urea-containing bioactive compounds. Whilst the conformational preferences of biaryl ureas have been extensively studied, very little attention has been paid to alkylated analogues. We carried out a systematic study of N-aryl (phenyl and pyridyl)-N′-cyclopentyl ureas with differing N-methylation patterns using Well Tempered Metadynamics at a semi-empirical level in implicit water (GBSA) using Well-Tempered Metadynamics to generate their conformational free-energy landscapes. Geometries and energetics of the most relevant configurations are further refined using DFT level of theory. Validation for the computation was obtained by synthesis of all 8 analogues followed by conformational studies by X-ray crystallography and NMR. These findings reveal that the methylation pattern significantly affects the conformational preference of the system. Most notably, N-phenyl-N′-cyclopentyl urea is shown to adopt both the trans–trans, and cis–trans conformations with equal energy and that the cis–trans conformation can be significantly stabilised by the presence of an internal hydrogen bond to the N′-hydrogen. This study will be of utility for the design of N-alkyl-N′-aryl ureas as drug candidates.

Conserved folding landscape of monomeric initiator caspases

The apoptotic caspase subfamily evolved into two subfamilies-monomeric initiators and dimeric effectors; both subfamilies share a conserved caspase-hemoglobinase fold with a protease domain containing a large subunit and a small subunit. Sequence variations in the conserved caspase-hemoglobinase fold resulted in changes in oligomerization, enzyme specificity, and regulation, making caspases an excellent model for examining the mechanisms of molecular evolution in fine-tuning structure, function, and allosteric regulation. We examined the urea-induced equilibrium folding/unfolding of two initiator caspases, monomeric caspase-8 and cFLIPL, over a broad pH range. Both proteins unfold by a three-state equilibrium mechanism that includes a partially folded intermediate. In addition, both proteins undergo a conserved pH-dependent conformational change that is controlled by an evolutionarily conserved mechanism. We show that the conformational free energy landscape of the caspase monomer is conserved in the monomeric and dimeric subfamilies. Molecular dynamics simulations in the presence or the absence of urea, coupled with limited trypsin proteolysis and mass spectrometry, show that the small subunit is unstable in the protomer and unfolds prior to the large subunit. In addition, the unfolding of helix 2 in the large subunit results in disruption of a conserved allosteric site. Because the small subunit forms the interface for dimerization, our results highlight an important driving force for the evolution of the dimeric caspase subfamily through stabilizing the small subunit.

Characterizing the Conformational Free-Energy Landscape of RNA Stem-Loops Using Single-Molecule Field-Effect Transistors.

We have developed and used single-molecule field-effect transistors (smFETs) to characterize the conformational free-energy landscape of RNA stem-loops. Stem-loops are one of the most common RNA structural motifs and serve as building blocks for the formation of complex RNA structures. Given their prevalence and integral role in RNA folding, the kinetics of stem-loop (un)folding has been extensively characterized using both experimental and computational approaches. Interestingly, these studies have reported vastly disparate timescales of (un)folding, which has been interpreted as evidence that (un)folding of even simple stem-loops occurs on a highly rugged conformational energy landscape. Because smFETs do not rely on fluorophore reporters of conformation or mechanical (un)folding forces, they provide a unique approach that has allowed us to directly monitor tens of thousands of (un)folding events of individual stem-loops at a 200 μs time resolution. Our results show that under our experimental conditions, stem-loops (un)fold over a 1-200 ms timescale during which they transition between ensembles of unfolded and folded conformations, the latter of which is composed of at least two sub-populations. The 1-200 ms timescale of (un)folding we observe here indicates that smFETs report on complete (un)folding trajectories in which unfolded conformations of the RNA spend long periods of time wandering the free-energy landscape before sampling one of several misfolded conformations or the natively folded conformation. Our findings highlight the extremely rugged landscape on which even the simplest RNA structural elements fold and demonstrate that smFETs are a unique and powerful approach for characterizing the conformational free-energy of RNA.

Friends in need: How chaperonins recognize and remodel proteins that require folding assistance.

Chaperonins are biological nanomachines that help newly translated proteins to fold by rescuing them from kinetically trapped misfolded states. Protein folding assistance by the chaperonin machinery is obligatory in vivo for a subset of proteins in the bacterial proteome. Chaperonins are large oligomeric complexes, with unusual seven fold symmetry (group I) or eight/nine fold symmetry (group II), that form double-ring constructs, enclosing a central cavity that serves as the folding chamber. Dramatic large-scale conformational changes, that take place during ATP-driven cycles, allow chaperonins to bind misfolded proteins, encapsulate them into the expanded cavity and release them back into the cellular environment, regardless of whether they are folded or not. The theory associated with the iterative annealing mechanism, which incorporated the conformational free energy landscape description of protein folding, quantitatively explains most, if not all, the available data. Misfolded conformations are associated with low energy minima in a rugged energy landscape. Random disruptions of these low energy conformations result in higher free energy, less folded, conformations that can stochastically partition into the native state. Two distinct mechanisms of annealing action have been described. Group I chaperonins (GroEL homologues in eubacteria and endosymbiotic organelles), recognize a large number of misfolded proteins non-specifically and operate through highly coordinated cooperative motions. By contrast, the less well understood group II chaperonins (CCT in Eukarya and thermosome/TF55 in Archaea), assist a selected set of substrate proteins. Sequential conformational changes within a CCT ring are observed, perhaps promoting domain-by-domain substrate folding. Chaperonins are implicated in bacterial infection, autoimmune disease, as well as protein aggregation and degradation diseases. Understanding the chaperonin mechanism and the specific proteins they rescue during the cell cycle is important not only for the fundamental aspect of protein folding in the cellular environment, but also for effective therapeutic strategies.

Entropic tuning regulates the function of biomolecular machines

Biomolecular machines execute the majority of essential biological processes in cells. Central to the function of these machines are thermally activated, spontaneous structural rearrangements that are governed by the corresponding conformational free-energy landscape. Because of their principal importance, fully characterizing the free-energy landscapes of biomolecular machines and determining how ligands modulate these landscapes in order to drive and regulate the functions of these machines remains a grand challenge in biology.

Boosting the conformational sampling by combining replica exchange with solute tempering and well-sliced metadynamics.

Methods that combine collective variable (CV) based enhanced sampling and global tempering approaches are used in speeding-up the conformational sampling and free energy calculation of large and soft systems with a plethora of energy minima. In this paper, a new method of this kind is proposed in which the well-sliced metadynamics approach (WSMTD) is united with replica exchange with solute tempering (REST2) method. WSMTD employs a divide-and-conquer strategy wherein high-dimensional slices of a free energy surface are independently sampled and combined. The method enables one to accomplish a controlled exploration of the CV-space with a restraining bias as in umbrella sampling, and enhance-sampling of one or more orthogonal CVs using a metadynamics like bias. The new hybrid method proposed here enables boosting the sampling of more slow degrees of freedom in WSMTD simulations, without the need to specify associated CVs, through a replica exchange scheme within the framework of REST2. The high-dimensional slices of the probability distributions of CVs computed from the united WSMTD and REST2 simulations are subsequently combined using the weighted histogram analysis method to obtain the free energy surface. We show that the new method proposed here is accurate, improves the conformational sampling, and achieves quick convergence in free energy estimates. We demonstrate this by computing the conformational free energy landscapes of solvated alanine tripeptide and Trp-cage mini protein in explicit water.

Subsecond Time-Resolved Mass Spectrometry in Dynamic Structural Biology

Life at the molecular level is a dynamic world, where the key players─proteins, oligonucleotides, lipids, and carbohydrates─are in a perpetual state of structural flux, shifting rapidly between local minima on their conformational free energy landscapes. The techniques of classical structural biology, X-ray crystallography, structural NMR, and cryo-electron microscopy (cryo-EM), while capable of extraordinary structural resolution, are innately ill-suited to characterize biomolecules in their dynamically active states. Subsecond time-resolved mass spectrometry (MS) provides a unique window into the dynamic world of biological macromolecules, offering the capacity to directly monitor biochemical processes and conformational shifts with a structural dimension provided by the electrospray charge-state distribution, ion mobility, covalent labeling, or hydrogen-deuterium exchange. Over the past two decades, this suite of techniques has provided important insights into the inherently dynamic processes that drive function and pathogenesis in biological macromolecules, including (mis)folding, complexation, aggregation, ligand binding, and enzyme catalysis, among others. This Review provides a comprehensive account of subsecond time-resolved MS and the advances it has enabled in dynamic structural biology, with an emphasis on insights into the dynamic drivers of protein function.

(Invited) Computational Simulations of Selective Interactions between ssDNA and SWCNTs

Structure-specific interactions between single-stranded DNA oligos and SWCNTs form the basis of important selective sorting methods yet remain poorly understood. Among these interactions is the differential affinity of (ATT)4 for one of the two (7,5) enantiomers. To understand the nature of the affinity difference, we have performed molecular dynamics (MD) simulations followed by free energy (replica exchange) calculations over the temperature range of 290 to 727 K. Comparisons between the conformational free energy landscapes of (ATT)4 structure on right-handed and left-handed (7,5) SWCNTs (respectively known as M and P in stereochemical convention) suggest that (ATT)4 has more contact area with M than with the P enantiomer. However, the corresponding end-to-end distances of (ATT)4 fall nearly in the same region on the landscape. Our standard MD simulations for (ATTT)3 interactions with SWCNTs also suggest that this DNA sequence covers (8,4) SWCNT more effectively than (7,6) and (8,3) SWCNTs, which is consistent with our recent dye-quenching experiments.Other experiments have revealed that when guanines in ssDNA oligos are covalently functionalized to SWCNT sidewalls, the extent of functionalization can depend on the position of the guanine within the oligo. Central positions within long oligos (e.g. T15GT15) lead to greater functionalization than end guanine positions (e.g. T30G). We have performed both standard MD and free energy calculations to investigate the source of this effect. For efficient guanine functionalization of SWCNTs, guanine needs to remain adsorbed on the SWCNT surface. Our computational observations suggest that the central guanine within T15GT15 oligo remains adsorbed on the (6,5) SWCNT surface even when the DNA interacts with other neighboring DNA strands. Under these conditions, the end guanine in T30G is destabilized and tends to desorb from (6,5) SWCNT surface due to inter-strand interactions. Our free energy calculations also suggest that oligos with a central guanine assume more elongated DNA conformations at the surface of a (6,5) SWCNT than do oligos with end guanines. We suggest that DNA elongation is associated with more stable adsorption on the SWCNT surface, supporting efficient covalent functionalization of central guanines. Figure 1

Molecular Dynamics Free Energy Simulations Reveal the Mechanism for the Antiviral Resistance of the M66I HIV-1 Capsid Mutation.

While drug resistance mutations can often be attributed to the loss of direct or solvent-mediated protein−ligand interactions in the drug-mutant complex, in this study we show that a resistance mutation for the picomolar HIV-1 capsid (CA)-targeting antiviral (GS-6207) is mainly due to the free energy cost of the drug-induced protein side chain reorganization in the mutant protein. Among several mutations, M66I causes the most suppression of the GS-6207 antiviral activity (up to ~84,000-fold), and only 83- and 68-fold reductions for PF74 and ZW-1261, respectively. To understand the molecular basis of this drug resistance, we conducted molecular dynamics free energy simulations to study the structures, energetics, and conformational free energy landscapes involved in the inhibitors binding at the interface of two CA monomers. To minimize the protein−ligand steric clash, the I66 side chain in the M66I−GS-6207 complex switches to a higher free energy conformation from the one adopted in the apo M66I. In contrast, the binding of GS-6207 to the wild-type CA does not lead to any significant M66 conformational change. Based on an analysis that decomposes the absolute binding free energy into contributions from two receptor conformational states, it appears that it is the free energy cost of side chain reorganization rather than the reduced protein−ligand interaction that is largely responsible for the drug resistance against GS-6207.

Ensemble simulations: folding, unfolding and misfolding of high‐efficiency transfer‐messenger RNA pseudoknot

The PK1 pseudoknot, a well-characterized super secondary structure, from Aquifex aeolicus is a transfer-messenger RNA (tmRNA) that participates in the ribosomal trans-translation rescue mechanism. Extensive all-atom molecular dynamics simulations were conducted across a worldwide distributed computing network to study the folding, unfolding, misfolding and conformational plasticity of the twenty-two nucleotides PK1 pseudoknot, yielded nearly 140 µs of cumulative sampling time. Complex K-means clustering was utilized to identify twenty-seven conformational microstates among forty starting structures spanning from the extended unfolded state to the native fold. Moreover, post equilibration data, established at 7.5 ns, provided detailed depictions the conformational free energy landscape that encompass the native states, numerous intermediate and misfolded states and intricate folding mechanism of the pseudoknot motif. The PK1 pseudoknot displays structural similarities to the class of frameshifting pseudoknots and suggests a parallel between the trans-translational mechanism and frameshifting mechanism of our previous study of the potato lead roll virus RNA pseudoknot.

Metamorphic proteins: the Janus proteins of structural biology.

The structural paradigm that the sequence of a protein encodes for a unique three-dimensional native fold does not acknowledge the intrinsic plasticity encapsulated in conformational free energy landscapes. Metamorphic proteins are a recently discovered class of biomolecules that illustrate this plasticity by folding into at least two distinct native state structures of comparable stability in the absence of ligands or cofactors to facilitate fold-switching. The expanding list of metamorphic proteins clearly shows that these proteins are not mere aberrations in protein evolution, but may have actually been a consequence of distinctive patterns in selection pressure such as those found in virus-host co-evolution. In this review, we describe the structure-function relationships observed in well-studied metamorphic protein systems, with specific focus on how functional residues are sequestered or exposed in the two folds of the protein. We also discuss the implications of metamorphosis for protein evolution and the efforts that are underway to predict metamorphic systems from sequence properties alone.

Charting the Complete Thermodynamic Landscape of Gas Adsorption for a Responsive Metal-Organic Framework.

New nanoporous materials have the ability to revolutionize adsorption and separation processes. In particular, materials with adaptive cavities have high selectivity and may display previously undiscovered phenomena, such as negative gas adsorption (NGA), in which gas is released from the framework upon an increase in pressure. Although the thermodynamic driving force behind this and many other counterintuitive adsorption phenomena have been thoroughly investigated in recent years, several experimental observations remain difficult to explain. This necessitates a comprehensive analysis of gas adsorption akin to the conformational free energy landscapes used to understand the function of proteins. We have constructed the complete thermodynamic landscape of methane adsorption on DUT-49. Traversing this complex landscape reproduces the experimentally observed structural transitions, temperature dependence, and the hysteresis between adsorption and desorption. The complete thermodynamic description presented here provides unparalleled insight into adsorption and provides a framework to understand other adsorbents that challenge our preconceptions.

Comprehensive characterization of oligosaccharide conformational ensembles with conformer classification by free-energy landscape via reproductive kernel Hilbert space.

Oligosaccharides play versatile roles in various biological systems but are difficult to characterize from a structural viewpoint due to their remarkable degrees of freedom in internal motion. Therefore, molecular dynamics simulations have been widely used to delineate the dynamic conformations of oligosaccharides. However, hardly any methods have thus far been available for the comprehensive characterization of simulation-derived conformational ensembles of oligosaccharides. In this research, we attempted to develop a non-linear multivariate analysis by employing a kernel method using two homologous high-mannose-type oligosaccharides composed of ten and eleven residues as model molecules. These oligosaccharides' conformers derived from simulations were mapped into reproductive kernel Hilbert space with a positive definite function in which all required non-redundant variables for describing the oligosaccharide conformations can be treated in a non-biased manner. By applying Gaussian mixture model clustering, the oligosaccharide conformers were successfully classified by different funnels in the free-energy landscape, enabling a systematic comparison of conformational ensembles of the homologous oligosaccharides. The results shed light on the contributions of intraresidue conformational factors such as the hydroxyl group orientation and/or ring puckering state to their global conformational dynamics. Our methodology will open opportunities to explore oligosaccharides' conformational spaces, and more generally, molecules with high degrees of motional freedom.

Metal–Organic Frameworks: On the Role of Enthalpic and Entropic Contributions to the Conformational Free Energy Landscape of MIL‐101(Cr) Secondary Building Units (Adv. Theory Simul. 12/2020)

A major cause of defects during the synthesis of metal–organic frameworks (MOFs) is the thermostructural behavior of the precursors used. In article number 2000092, Vassiliki-Alexandra Glezakou, Matteo Salvalaglio, and co-workers use enhanced sampling methods to evaluate the flexibility of MIL-101(Cr) half-secondary building units (SBUs) in solution. They find that the interplay between enthalpy and entropy determines the most probable conformational state for such precursors, highlighting the importance of explicitly considering entropic effects when estimating the relative stability of different conformers of SBUs.

On the Role of Enthalpic and Entropic Contributions to the Conformational Free Energy Landscape of MIL‐101(Cr) Secondary Building Units

AbstractThe thermo‐structural behavior of metal‐organic framework (MOF) precursors is responsible for regulating the introduction of defects in MOF structures during synthesis. In this paper, factors affecting the flexibility of MIL‐101(Cr) half‐secondary building units (half‐SBUs) are evaluated in solution using enhanced sampling methods. In particular, entropic and enthalpic contributions to the conformational free energy landscape of isolated MIL‐101(Cr) half‐SBUs are calculated in water, in the presence and absence of ionic species (Na+ and F−), and in N, N‐dimethylformamide (DMF). This analysis leads to the observation that the interplay between enthalpy and entropy determines the most probable conformational state for half‐SBUs. This observation extends to the most relevant SBU intermediate, in which conformational entropy plays a key role in stabilizing configurations that differ from those found in the MIL‐101(Cr) crystal structure. The findings highlight the importance of explicitly considering entropic effects, associated with finite‐temperature sampling when estimating the relative stability of different conformers of SBUs.

D-Serine Potently Drives Ligand-Binding Domain Closure in the Ionotropic Glutamate Receptor GluD2

Despite their classification as ionotropic glutamate receptors, GluD receptors are not functional ligand-gated ion channels and do not bind glutamate. GluD2 receptors bind D-serine and coordinate transsynaptic complexes that regulate synaptic plasticity. Instead of opening the ion channel pore, mechanical tension produced from closure of GluD2 ligand-binding domains (LBDs) drives conformational rearrangements for non-ionotropic signaling. We report computed conformational free energy landscapes for the GluD2 LBD in apo and D-serine-bound states. Unexpectedly, the conformational free energy associated with GluD2 LBD closure upon D-serine binding is greater than that for AMPA, NMDA, and kainate receptor LBDs upon agonist binding. This excludes insufficient force generation as an explanation for lack of ion channel activity in GluD2 receptors and suggests that non-ionotropic conformational rearrangements do more work than pore opening. We also report free energy landscapes for GluD2 LBD harboring a neurodegenerative mutation and demonstrate selective stabilization of closed conformations in the apo state.